Vortices Freeze like Window Glass: the Vortex Molasses Scenario 
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We overview several recent experimental and numerical observations, which are at odds with the 
Vortex Glass theory of the freezing of disordered vortex matter. To reinvestigate the issue, we 
performed numerical simulations of the overdamped London - Langevin model, and use finite size 
scaling to analyze the data. Upon approaching the transition the initial Vortex Glass type criticality 
is arrested at some crossover temperature. Below this temperature the timescales continue growing 
very quickly, consistent with the Vogel-Fulcher form, while the spatial correlation length ^ stops 
exhibiting any observable divergence. We call this mode of freezing the Vortex Molasses scenario. 



PACS numbers: 74.60. Ge, 74.25. Dw 

The influence of disorder on vortex matter is one of the 
most paradigmatic problems. The vortex lattice, formed 
in clean systems, is inherently unstable towards a less or- 
dered state even for infinitesimally small disorder ^ . At 
small magnetic fields, or equivalently, at weak disorder, 
a dislocation free phase is emerges, which thus retains 
a topological order. Small angle neutron scattering [gj, 
and Bitter decoration experiments H] seem to support 
this picture. The theoretical foundation for such a Bragg 
Glass was provided by scaling arguments H, and varia- 
tional calculations |q| . Its cornerstone is the logarithmic 
behaviour of vortex correlations at large distances. Nu- 
merical simulations also reported a strongly suppressed 
dislocation density p] , and confirmed the logarithmic be- 
haviour of correlations M below a critical field strength. 
It is noteworthy, however, that the largest scale imaging 
studies S, did not find evidence for logarithmic correla- 
tions, thus the details of the dislocation free regime are 
still subject to discussion. 

We also understand the influence of increasing fields, 
or disorder. The key phenomenon here is the appearance 
of dislocation loops, accompanied by the entanglement of 
vortices ^,|l^ . Experimental support for this idea is the 
sharp enhancement of the critical current from magne- 
tization measurements |ll|-|l4|, the rapid destruction of 
the Bragg peaks in neutron scattering ^, and the pro- 
nounced dips in the electric field - current density, or E-J 
curves llq] . Numerical studies found evidence for disloca- 
tion loops destroying a quasi-ordered state in frustrated 
XY models [|6| , in Lawrence-Doniach representations [^ , 
and in realistic London-Langevin approaches R]. 

The nature of the high field phase is still very much 
in debate. The thermally assisted flux flow (TAFF) pic- 
ture predicts that vortices move in bundles, and over- 
come barriers via thermal excitations. This destroys su- 
perconductivity because the linear resistivity assumes a 
finite value, governed by an activated temperature de- 
pendence R{T) r^ Ro exp(-C//T) Q. An influential al- 
ternative was put forward in the form of the Vortex Glass 
(VG) theory |18 19]. The proposed Vortex Glass phase 



is distinguished by an unbounded distribution of barrier 
heights. This results in the vanishing of the linear resis- 
tivity, thus restoring superconductivity, and inherently 
non-linear E-J characteristics. Numerical support for this 
picture emerged from the study of the isotropic Gauge 
Glass model, ignoring the effects of screening [g^ . Exper- 
imental confirmation soon followed, on heavily twinned 
YBCO films [Q. A key evidence was provided by ob- 
serving the scaling of a crossover current Jx [21) . 

Recent experimental and numerical work, however, has 
raised new questions about the Vortex Glass picture. 

i) The values of the correlation length, creep and dy- 
namical exponents, i', /i and z, respectively, seem to de- 
pend on temperature, current, and sample quality ^^ in 
a very nonuniversal way. ly was found between 1.3 — 2, z 
between 3.1 — 6.5, and fj, between 0.2 — 0.5. 

m) The above values of v and z are much higher than 
their mean field values, indicating that the lower critical 
dimension might be close to 3. It is already accepted that 
there is no finite temperature VG phase in 2D p^ . 

Hi) Recent experiments in completely untwinned 
YBCO samples found that the E — J curves remained 
completely linear down to the lowest measurable values 
of the current. Correspondingly no scaling behaviour of 
the E — J curves were found |^,^ . This suggests that 
in previous works the twin boundaries might have played 
the role of extended defects, and in fact the observed 
scaling behaviour was that of the Bose Glass. 

iv) When the twin boundaries were removed in YBCO 
samples, the crucial crossover current J^ was found to 
saturate, instead of exhibiting a scaling behaviour ||26|| . 

v) Recent numerical works reported that when a finite 
London screening length A was restored into the previ- 
ously studied Gauge Glass models, the finite tempera- 
ture Vortex Glass transition disappeared p^-[29|]. Now, 
close to the transition the vortex correlation length is 
supposed to diverge, thus exceeding A. Therefore the ul- 
timate transition region is always in this finite A regime. 

We conclude that there is an emerging body of evi- 
dence, which is inconsistent with the Vortex Glass pic- 



ture. Motivated by this inconsistency, in this paper we 
explore analogies to an other widely studied glass transi- 
tion, and investigate the possibility that Vortices freeze 
like the window glass: the Vortex Molasses scenario. 

We realize that there is not a uniquely accepted theory 
of the window glass transition. Therefore we construct 
the Vortex Molasses (VM) scenario only from those ele- 
ments, which are common among the different theories: 
i) a very rapid freezing of the dynamics, with diverg- 
ing timescales, characterized by the Vogel - Fulcher law: 
r ~ exp[l/(T — Tg)]; ii) the possible divergence of the 
spatial correlation length ^ is rendered unobservable by 
this rapid freezing. 

We note that among the early alternative propositions, 
some emphasize the entanglement of vortices in the pres- 
ence of disorder, in analogy to polymer glasses [ pO[. Also, 
extensions of the TAFF theory were constructed |3l| . Fi- 
nally, a Vortex Slush picture has been proposed, viewing 
the glass of vortices as a viscous liquid, driven by the 
remnants of the first order melting transition ||32| . 

We start by overviewing ways to distinguish between 
the VG and VM scenarios. First, the predicted tempera- 
ture dependence of the resistivity differs: in the VG the- 
ory the resistivity vanishes as p{T) — pvg\T — Tg\''''^~^\ 
whereas in the VM scenario one expects the resistivity 
to follow the Vogel-Fulcher law p{T) = pvMexp[—l/{T — 
Tg)] 133]. However it is hard to achieve decisive distinc- 
tion between these forms, as the resistivity exponent in 
the VG theory is large: u{z - 1) r^ 5 - 7 ^. 

The second method is more promising. Equating the 
current related free energy of a correlated volume to the 
thermal energy yields the above mentioned crossover cur- 
rent density scale J^ = cfcsr/($oC^) |19|. Above the 
transition temperature the low current linear E ^ J oi 
the viscous liquid is expected to cross over at J^ to the 
Bardeen - Stephen form at high currents. In the VG the- 
ory the correlation length ^ diverges: £, ^ {T — Tg)^'^ ■ 
Correspondingly the crossover current scale Jx collapses 
asJxiT) - [T-TgY". In contrast, in the VM picture the 
criticality of the correlation length is unobservable^ thus 
Jx does not collapse. As mentioned, while in twinned 
samples Jx collapsed, in untwinncd YBCO Jx saturated 
at some finite value upon approaching the transition p6[ . 

In this paper we report numerical simulations of a real- 
istic, London-type model for driven vortices, governed by 
overdamped dynamics in order to distinguish between the 
above two scenarios. Previous studies on the gauge glass 
already indicated a breakdown of the VG picture p7|-p9[ , 
however that model is rather simplified. For instance it 
is isotropic, whereas in real vortex matter the external 
field definitely introduces a strong anisotropy. Thus it 
remains an open question, whether the gauge glass ade- 
quately describes the vortex matter, making our realistic 
simulations necessary. The Langevin equation describing 
the overdamped motion of vortices reads 




FIG. 1. Tlie current dependence of tlie resistivity 
p{T,I)/pBS at T/Tc = 0.74. Inset: The temperature depen- 
dence of the low - current ohmic resistance p{T, J— > 0)/pbs, 
for L = 8 systems. Dark gray: fit to the Vogel - Fulcher 
form: p{T) = poexp[-ro/(r - Tg)], light grey: VG fit: 
p{T) = po[(T-Tg)/Tg]^-. 
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where p labels the vortices with coordinates R^, Ciz,t) 
is the Langevin noise, and F^ is the Lorcntz force. The 
Hamiltonian H is constructed on the basis of the Lon- 
don theory. Its derivative decomposes into three forces: 
the pairwise interactions, the single- vortex bending force, 
and the pinning. For more details of the method, see Ref. 
(Q. The trustworthiness of our code was demonstrated 
by the quantitatively correct reproduction of the phase 
diagram of disordered YBCO H . In our model we do not 
take into account vortex-loops whose effect on melting is 
still controversial 34 1. 

There are several length scales in the model. To avoid 
observing some crossover instead of the asymptotic be- 
haviour, we chose the characteristic microscopic length 
scales small and close to each other. We use A/^ — 4, 
a large magnetic field of H/Hc2 = 0.2, to make oq sa A, 
and finally we made the system isotropic by choosing 
e = 1. The 100-500 vortex elements produce good self- 
averaging, so a reasonable statistics was achieved by aver- 
aging over 10-20 disorder realizations. Simulated anneal- 
ing was employed to generate the starting configurations. 

In Fig. 1 we show the typical behaviour of the differen- 
tial resistivity p{T,I) = dE/dJ, normalized to pss, the 
Bardeen-Stephen value. At high currents p{T, J) = pss 
as it should. After a pronounced drop with decreasing 
current, p{T, J) fiattens at low J, clearly indicating an 
ohmic behaviour: we are in the Vortex Liquid regime. 

In our temperature sweeps p{T) drops by two orders 
of magnitude upon approaching the freezing transition. 
The inset of Fig. 1 shows its temperature dependence. We 
also exhibit a Vogel - Fulcher fit, and a power law fit. 
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FIG. 2. The temperature dependence of R{T) for systems 
with L = 4, 6, 8. 



As expected, both fits are comparable, and thus do not 
distinguish between the VM and VG theories. 

However, useful information can be extracted from the 
R(T) runs via finite size scaling. In Fig. 3 we show the 
results for 4'^, 6'^, and 8'^ systems. Upon cooling, the 
system exhibits increasing finite size sensitivity down to 
T/Tc ~ 0.70, which typically indicates the approaching 
of a phase transition. Remarkably, however, below this 
temperature range the finite size sensitivity decreases, as 
if the criticality is arrested. 

To view this from a different perspective, we follow 
Young et al. |g^ by studying an analogue of the Binder 
- ratio, log[p{L) / p{L')] / log[L / L']. Far from a transition, 
where finite size sensitivity is small, this should be close 
to 0. Approaching a transition the increased finite size 
sensitivity is signalled by the data on different system 
sizes splaying out. Eventually, however, they come to- 
gether and cross atT = Tq. In Fig. 4 we show this ratio 
for our model. As the temperature decreases, the ini- 
tial splaying shows an impending transition. The curves, 
however, do not cross. Instead, they turn back up: the 
transition is arrested. This again signals the decrease of 
finite size sensitivity, which is most readily interpreted 
as the initial increase of the correlation length ^ being 
arrested around T/Ta ~ 0.70. 

The best measure of i^ is via the crossover current Jx , 

— 1/2 

plotted in Fig. 5. As ^ ^ J^; , decreasing Jx indicates 
increasing correlation length. However Jx, and thus ^, 
saturates with decreasing T, around T/Tc '^ 0.68, in 
quantitative accord with the finite size scaling. Note that 
a very similar flattening of Jx was observed in untwinned 
YBCO pq]. All these three tests can be interpreted as 
follows. Upon decreasing the temperature a Vortex Glass 
criticality starts to develop. However this critical be- 
haviour gets arrested around T/Tc — 0.69 ± 0.01, and 
crosses over to a Vortex Molasses criticality. This is char- 
acterized by a rapid, Vogel - Fulcher type decrease of the 




FIG. 3. The logarithmic resistance ratio as a function of 
the teMpi e|raturc. 




FIG. 4. The temperature dependence of the crossover cur- 
rent Ix, for L = 4,6, 8. 



resistivity, but at the same time an essentially noncritical 
behaviour of the correlation length ^. 

The above results established that the freezing tran- 
sition is unlikely to be governed by Vortex Glass the- 
ory, but rather it looks more like a window glass tran- 
sition. However there isn't a single theory agreed upon 
by the window glass community. For a review of dif- 
ferent approaches, see Ref. [Q. Some theories propose 
that the correlation length diverges as a power law, but 
the freezing of the dynamics is so rapid, that it ren- 
ders this divergence unobservable. Others believe that 
in fact ^ does not diverge at all, it remains noncritical 
even on the longest time scales. Finally there are theo- 
ries which envision that there is no true transition at any 
finite temperatures, but a rapid, continuous increase of 
the viscosity, diverging only at T = 0. This latter view 
was imported to the vortex problem by Ref. [2^. Set- 
ting up the scaling theory for finite size systems with a 
correlation length ^ ~ T~'^ gives for the nonlinear E-J 
relation: E/{JR) = E{J/T^+^'', L^/^T) , where £: is a 



universal scaling function. Adopting the accepted defini- 
tion of the crossover current density, E/{JxR) = 2 yields: 
Jx = T^+'^''f{L^/"T) , where / is an other universal 
function. Ref. ||2^ finds Jx to be a universal function 
of L^/'^T, with u = \. To test this proposition, we also 
plotted JxlT^^"^" as a function of L^/T. However we 
found no universal dependence whatsoever. This clearly 
eliminates the possibility of a T = fixed point govern- 
ing the freezing behaviour of the London model. Thus, 
remarkably, regarding the freezing transition the gauge 
glass and the realistic vortex simulations give qualita- 
tively different results. 

We now understand that structural (window) glasses 
and systems with quenched disorder often behave quite 
similarly |3q ]. Their glassy phase exhibits different ag- 
ing phenomena [|37| . Measuring the two time correlation 
functions ^^ of vortices, and comparing to the predicted 
power law relaxation forms would be a constructive test 
of the Vortex Molasses. 

A word on the appropriateness of models with strong 
screening. Ref. [g^ recalls that in the analogous 3d XY 
model hoth the screening length A and correlation length 
^ diverge, when T^ is approached from below. In the criti- 
cal region the exponent of A is half of ^'s, and hence close 
enough to Tc the proper characterization of the system 
should involve strong screening. Above Tc the screening 
length on macroscopic scales is infinite, as we are in the 
Vortex Liquid. On the scale of intervortex separation, 
however, A equals its bare value. In general, the pres- 
ence of the other vortices generates a renormalized, scale 
dependent X{x). Whether the model is in the strong or 
weak screening limit, will then be determined by A(^) 
being greater or smaller than ^. By invoking that the 
critical behaviour around Tc is typically symmetric, and 
that below Tc we are in the strong screening limit, we 
expect A(^) < ^, i.e. the screening remaining essential 
for understanding the physics of the model, the starting 
point of our simulations. 

In conclusion, we collected several numerical and ex- 
perimental results, which are at odds with the Vor- 
tex Glass theory of the freezing of the disordered vor- 
tex matter. Previous confirmations of the VG the- 
ory were reinterpreted in terms of twin boundaries and 
proper accounts of the screening. To reinvestigate the is- 
sue, we performed careful numerical simulations of the 
overdamped London - Langevin model, and used fi- 
nite size scaling to analyze the data. We found that 
upon approaching the transition the initial Vortex Glass 
type criticality is arrested at some crossover temper- 
ature, where the vortex correlation length catches up 
with the screening length. Below this temperature the 
timescale continues growing very quickly, consistent with 
the Vogel-Fulcher form, while the spatial correlation 
length ^ stops exhibiting any observable divergence. We 
call this mode of freezing the Vortex Molasses scenario. 
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